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Constraints on the whole spectrum of lepton flavor violating vertices are shown in the context of the 
standard two Higgs doublet model. The vertex involving the e — r mixing is much more constrained than 
the others, and the decays proportional to such vertex are usually very supressed. On the other hand, 
bounds on the quark sector are obtained from leptonic decays of the B® s mesons and from AM B o . We 

emphasize that although the B% — B d mixing restricts severely the d — b mixing vertex, the upper bound 
for this vertex could still give a sizeable contribution to the decay B® — > j/p respect to the standard 

l/~) , model contribution, from which we see that such vertex could still play a role in the phenomenology. 
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1 Introduction 

O 

Many extensions of the standard model leads naturally to flavor changing neutral currents (FCNC) in quark 
and lepton sectors. This is the case of models with an extended Higgs sector. However, owing to the high 
JL , suppression imposed by experiments, several mechanisms has been used to get rid of them, such as discrete 

symmetries PP, permutation symmetries PJ, and different textures of Yukawa couplings |3J. Notwithstanding, 
the increasing evidence on neutrino oscillations seems to show the existence of mass terms for neutrinos as 
well as of family lepton flavor violation (LFV) [3] . Such fact has inspired the study of many scenarios that 
predict LFV processes as in the case of SUSY theories with R-parity broken jS], SU(5) SUSY models with 
right-handed neutrinos [H], models with heavy Majorana neutrinos [7], and multi-Higgs doublet models with 
right-handed neutrinos for each lepton generation ;$;. On the other hand, LFV in the charged sector has 
been also examined in models such as SUSY GUT |j]|, and the two Higgs doublet model (2HDM) |HllTT)j. 

In the charged lepton sector, searches for FCNC have been carried out through leptonic and semilep- 
tonic decays of K and B mesons ^J, as well as purely leptonic processes |12|. On the other hand, some 
collaborations plan to improve current upper limits of some LFV decays by several orders of magnitude, 
by increasing the statistics |1H| . Other possible sources of improvement lies on the Fermilab Tevatron and 
LHC by means of LFV Higgs boson decays. Further potential sources to look for Higgs mediated FCNC lie 
on the muon colliders. It is because they have the potentiality to produce Higgs bosons in the s— channel, 
with substantial rate production at the Higgs mass resonance ^3] • From the theoretical point of view, since 
Higgs Yukawa couplings are usually proportional to the lepton mass, they give an important enhancement 
to cross sections with Higgs mediated s— channels, respect to the ones in an e + e~ collider. 

As for the quark sector, it is well known that the data from K° — K and B° — B mixing put severe 
bounds on the flavor changing couplings involving the first family ^] ^] . Indeed, this fact was one of the 
motivations to implement a discrete symmetry in the 2HDM in order to suppress FCNC effects 1 . This fact 
in turn motivated the construction of a parameterization in the 2HDM in which the FC vertices involving 
the first family are neglected, and assume that the only non- vanishing couplings are \u, \b |17) . Based on 
this assumption, constraints on A^t and A^j, from B° — B and lower bounds on m^± from the CLEO data of 
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b — > 57 have been estimated [TB| . In this scenario, the tree- level diagrams for B°, — B d are neglected and the 
box diagrams involving one charged Higgs boson in the loop become dominant, leaving Xtt, Xbb and to#± as 
the only free parameters in the process. However, although Ref. ^j] found that relatively light charged 
Higgs bosons are still allowed, they also found that very heavy charged Higgs bosons are still permitted and 
even required if there is a significant relative phase between X tt and Xbb- Wherever a very large value of 
m H ± is allowed, it opens the possibility of having dominant or at least competitive tree level diagrams even 
with highly suppressed values of the couplings involving the first family, especially in the case in which at 
least one of the neutral Higgs bosons is kept light. Inspired in this idea, we shall assume in this paper that 
the tree level diagram is dominant. 

On the other hand, in a recent previous work |1(J|. some constraints on LFV have been found in the 
framework of the two Higgs doublet model with flavor changing neutral currents. Specifically, bounds 
on the vertices £^-r,£eT, £^,£tt, were obtained based on the g — 2 muon factor and the leptonic decays 
/i — > ej, t — > fififi, t — > fj.'j. Additionally, upper limits on the decays r — ► e"f and r — » eee were estimated, 
finding them to be highly suppressed respect to the present experimental sensitivity. The purpose of this 
work is on one hand to complete the information about the spectrum of the LFV matrix in the lepton sector, 
and on the other hand restrict some vertices involving the first family of the quark sector, and show that 
such vertices could still play a significant role in the phenomenology. Combining bounds on the quark and 
lepton sector we can predict upper bounds for leptonic decays of the B° mesons. 

2 Constraints in the lepton sector 

We shall work in the context of the two Higgs doublet model (2HDM) with flavor changing neutral currents, 
the so called model type III. We shall neglect possible relative phases between the FC vertices. The leptonic 
Yukawa couplings read 
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where H°(h°) denote the heaviest (lightest) neutral CP— even scalar, and ^4° is a CP— odd scalar. E refers to 
the three charged leptons E = (e, /i, r) and Me, £,e are the mass matrix and the LFV matrix respectively, 
a is the mixing angle in the CP— even sector. We use the parameterization in which one of the vacuum 
expectation values vanishes. 

The decays needed to obtain our bounds are given by 
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Observe that the decays containing two identical particles in the final state possess interferences involving 
the pseudoscalar Higgs boson, while the decays with no identical leptons in the final state do not con- 
tain interference terms with the pseudoscalar. On the other hand, in the calculation of the decay width 
r (t~ — » /i + /i~e~), we neglect diagrams containing the vertex £ er and keep only the ones proportional to 
fp r , we make this approximation because a previous phenomcnological analysis shows a strong hierarchy 
between these mixing vertices ^U| (|£er| << |£mt| by at least five orders of magnitude). 
The corresponding experimental upper limits for these rare processes are |18| 
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2.1 Bounds on £ Me and £ ee 

In a previous work |1U| . the LFV vertices coming from the 2HDM type III, were constrained by using several 
pure leptonic processes, the following bounds for the LFV vertices were found 
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Such constraints are valid in most of the region of parameters. Since we intend to complete the analysis 
made in |1U| . we shall make the same assumptions which we summarize here for completeness. We settle 
to/jO sa 115 GeV, and m^a > m h o. In order to cover a very wide region of parameters, we examine five cases 
for the remaining free parameters of the model |10j 

1. When m H o ~ 115 GeV. 

2. When m H a ~ 300 GeV and a = ir/2. 

3. When rxijja is very large and a — 7r/2. 

4. When m H o ~ 300 GeV and a = tt/4. 

5. When m^o is very large and a = it /A. 



and 



For all those cases the value of the pseudoscalar mass is swept in the range of m^a > 115 GeV. 

The vertex £ 2 can be constrained by combining the existing limits on £ 2 given in Eqs. |J3J, 
experimental upper limit on the branching ratio Br(r~ —* ^t~/i~e + ) given by Eq. J2J- Alternatively, we 
can constrain the same vertex from the decay r~ — v /i + /i~e~. The upper limits on £ 2 obtained from both 
decays are illustrated in table (JTJ for the five cases explained above. We should observe that the upper limits 



case 


from t —> /! /! e + 


from r — ► fi + fi e 


1 


Q e < 5.59 x 10- a 


Qe < 1-0 X 10- 2 


2 


e, e < 1.5 x 10- 1 


£ 2 e < 2.7 x 10- 1 


3 


unconstrained 


unconstrained 


4 


e ue < 1.35 x 10- 2 


e ue < 2.43 x 10- 2 


5 


Q e < 1.67 x 10- 2 


g e < 3.0 x 10- 2 



Table 1: Bounds on the mixing vertex £ 2 , based on the processes r 
the five cases cited in the text. 
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are less restrictive than the ones coming from r — > ^i /i e + . However, both 
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sets of constraints lie roughly on the same order of magnitude. From tabled we can extract a quite general 
bound for the vertex & ie 

& < 1.5 x lCr 1 , (4) 

valid for most of the region of parameters 1 . It is worth saying that other restrictions on this vertex can be 
gotten from u. — > ej or r — > ej assuming that only the diagrams with a muon in the loop contribute, instead 
of the tau as customary. However, bounds obtained this way are much less restrictive. 
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Figure 1: Contourplots for the five cases cited in the text in the £ ee — m^o plane, based on the process 
r~ — > /i~e + e~. On left: Case 1 (dotted line), case 4 (dashed line) and case 5 (solid line). On right: Case 2 
(solid line) and case 3 (dashed line). 



On the other hand, we can get contraints on the vertex £ ee by combining the already mentioned bounds 
on ^ 2 and the upper experimental constraints for the decay Y (t~ — ► /i~e + e~) of Eq. (J2J. Since the factor 
£ ee cannot be factorized in contrast to the case of £ Me , we extract its bounds in the form of contourplots 
in the £ ee — m^o plane, as shown in Fig. ft}. Additionally, we write in table (J2J, the constraints obtained 
for m^o very heavy and for m A a ss 115 GeV. From table © we can extract general constraints for £ ee , the 



Case £ ee (m^o very heavy) £ ee (m^o ~ 115 GeV) 


1 < 9.75x10"^ <6.89xl(T 2 


2 < 5.1X10" 1 < 7.41 x IQ-' 1 


3 unconstrained unconstrained 


4 < 1.5 x 1CT 1 <7.54xl0- 2 


5 ^UxlO" 1 <7.53xl(T 2 



Table 2: Bounds for the mixing matrix element £ e e, for m^o ~ 115 GeV and for m^o very heavy. Such 
constraints are based on the bounds on £ MT and the upper limit for the decay width Y (t~ — ► fi~e + e~). 



general bounds read 



\U\ <5.1x 10- 1 ; l^eel <7.54x 10" 



for tUao ~ 115 GeV and for tua very heavy respectively. We emphasize again that this prediction is valid 
in most of the region of parameters but fails in the case 3 cited above, i.e. when m H o is very large and 
a = it/2. 

Finally, we make a prediction about the upper limit for the branching ratio of the process r~ — » /i + e - e~, 
based on the limits on £ efi shown in table^^d the limits on £ eT shown in Eqs. (|3J), the results are collected 
in table |31 We see that the upper limits shown in table © are at least ten orders of magnitude smaller than 
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the present experimental upper limit Br (t~ — ► /i + e~e~) < 1.5 x 10~ 6 GeV, (except for the third case). 
In addition, table Q] collects the results for the general upper limits of three leptonic decays involving the 
vertex £ er . The strong suppression of these processes might be anticipated because of its proportionality to 
£e T which is much more restricted than the others |10j . 



Case 


Br(r — > fi + e e ) 


1 


< 9.5 x l(T iw 


2 


< 3.2 x 1(T 17 


3 


Unconstrained 


4 


< 1.3 x 1(T 17 


5 


< 1.4 x 10~ iY 



Table 3: Upper limits for the branching ratio Br (r — > e e /i + ), based on the contraints obtained for the 
LFV vertices £ Me and £ eT . The experimental upper limit is 1.5 x 10~ 6 



3 Constraints in the quark sector 

We shall obtain constraints on the quark sector by using the experimental information from B® lepton 
decays and Am B o. The B® measurements are dominated by the asymmetric B factories [T^J, while the B® 
measurements come from hadron colliders [201 • At the tree level, the decays B® — > IV depend on the product 
€?l J £db an< ^ the pseudoscalar Higgs boson mass only. In the framework of the 2HDM-III they are 
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Where /g represent the B^ meson decay constant whose value has been taken from |5J. The present 
experimental upper bounds for these decays at 90% C.L. are ^B] 
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Consequently, we can get upper limits for this products of mixing vertices by using the upper bound for 
these decays. In particular, from B d — * fj/r we find 

&•& £ (5-45 x 10- 15 GeU- 4 ) rn^o . 



Predictions 


Experim. limits 


Br(r -► e~7) < 6.6 x 10~ iB 


2.7 x 10 _e 


Br (t~ -> e+e-e") < 2.2 x 10~ 17 


2.9 x 10~ B 


Br (r -> ^+e~e-) < 3.2 x 10~ 1Y 


1.5 x 10~ e 



Table 4: Upper limits predicted for some lepton decays. All of them are highly suppressed respect to the 
current experimental upper limit. 
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On the other hand, since we have a lower bound on the mixing vertex £f; T we can obtain an upper bound 
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lb alone 
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We can obtain a similar bound for the product £ei-£db based on the upper limit for Bd — > er. Nevertheless, 
a better bound for this product is obtained by combining Eqs. ©,10. The same situation occurs for the 
product ^eu^db' sm ce the present bounds on the decays Bd — ► e/i cannot provide a better constraint for this 
product that the bound obtained from Eqs. ©, ©■ 

Furthermore, from the upper limit for Br (B® — ► e/i) we find an upper bound for the product £,^ b C, 



sb^e^j. 



e sb Z% < 2.38 x 10- I7 m\ 



As we see from Eq. ||Sj and the bounds in this section, the latter blow up rapidly when m A o grows. Indeed, 
for heavy values of the pseudoscalar mass, the one loop contributions could be sizeable [22| introducing more 
free parameters to the model. Nevertheless, if we consider a quite light pseudoscalar i.e. a mass no far from 
the electroweak scale, the tree level contribution is dominant and the bounds above are quite restrictive. 
Table shows some typical values for the upper limits above for 115 GeV< m A o < 250 GeV. 



TO A (GeV) 115 200 250 


<^ 6 (xl0- 3 ) 1.25 11.44 27.92 


&&■(>< 10- 1Y ) 3.46 31.69 77.35 


^ M (xl0- 4 ) 1.87 17.16 41.88 


€ibSw( xl0 ~ a ) °- 42 3 - 82 9 - 32 



Table 5: Bounds for £| b , £^ 6 and their products with leptonic vertices for pseudoscalar bosons lying roughly 
in the electroweak scale. 



3.1 Constraints from B% — B d mixing 

Now let us use the AM B o parameter to constrain the vertex £db involving the first family. The invariant 

d 

amplitude at the tree level for B q — B® mixing in our model is given by 



(B° q \H w \B° q ) = --^R 2 qm (B° q \bqbq\B* q ) - -^R 2 qbh (B q \bqbq\B* q ) 



'H 
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where R q bh are the coefficients of the Feynman rules with q — d, s. In terms of the operators 0^ F=2 defined 
in ^J E] > we have that 



(B"\0^\W) = B B (W\0^\W 



'VIA, 



(8) 



where VIA denotes the Vacuum Insertion Approximation and Bb is the vacuum saturation coefficient. The 
operators that we need in our case are 



Op = (bq)(bq) ; Op' = (&75«)(&75<z)- 



(9) 



In addition, based on the expressions shown in section VI of the first of Refs. ^U, we introduce the factors 
Mg and Mp in terms of the AF = 2 matrix elements of the only two operators which do not vanish in the 
vacuum 
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Case \£ db \ (m A = 115 GcV ) |£ di ,|( m A very large ) 


1 2.69 x 1(T B 7.37 x 1(T 6 


2 2.86 x 1(T H 19.22 x 10~ B 


3 2.89 x 10~ B unconstrained 


4 2.77 x 1(T B 9.73 x 10~ B 


5 2.79 x 10~ B 10.42 x lO" 8 



Table 6: Constraints on |^fc| for rriA = 115 GeV and tua very large, based on the data of AM Bd 
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Using Eqs. (f7|)- (|llf) . we find that the contribution from new physics to the parameter AM B o reads 
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We shall estimate bounds for £<j& by using AM B o coming from 2?^ — B^ mixing. The predictions for 
the standard model (SM) AMsm and the experimental value AMexp have been taken from ^H] by using 
symmetrical uncertainties for the sake of simplicity 



AM SM 
AM E xp 
The maximum room for the new physics reads 
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(14) 



where AM EXP and AMg M represent the central values of AMexp and AMsm respectively. Furthermore, 
Esm and Eexp are the uncertainties associated to the standard model and experimental estimations re- 
spectively. All of them are given by Eq. H13|) . So we can constrain £<$ based on the values of AMsm and 
AMexp for the B® meson. We found the results shown in table [5] for the five cases, when we consider 
ttia = 115 GeV and m^o very large. 

Now, using the bounds on Q b obtained from the AM B o and combining them with the allowed values 
for the vertex £ M/i in Eq. JSJ) , we shall predict the maximum contribution of this new physics to the decay 
B® — > jlp,. Such decay has already been considered in the literature in the framework of the two Higss 
doublet model with and without FCNC j23|. The predicted upper bounds for this decay for the five cases 
explained in the text, are displayed in table Hand the SM prediction which was calculated by avoiding the 
big uncertainties of f B o is given by [23 

Br (B d -» nfi) gM = 1 x 10- 10 

Taking into account the upper limits of \S,db\ we obtain that the tree level contribution to this process 
coming from the 2HDM can be comparable and even dominant respect to the SM contribution. It shows 
that, although the mixing vertices involving the first generation are highly suppressed, is still possible for 
them to play a role in the phenomenology. 
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Case Br(Bd — > /J,fi) (tua = 115 GcV ) Br(Bd — ► pt/x) ( to^ very large ) 


1 2.14 xlO~ 8 1 x 10~ 8 


2 2.3 x 10~ s 1 x 1CT* 


3 2.32 x 1CP S unconstrained 


4 2.21 x 1CT 8 1 x 10~ 8 


5 2.23 x 1CT 8 1 x 1CT 8 



Table 7: Upper limits for the branching ratio Br (B e i 
allowed values for £ w . 



/i/i) based on the upper limits for £db and the 



4 Conclusions 

We have found constraints on the whole spectrum of the mixing matrix of leptons, by using purely leptonic 
processes. A strong hierarchy between the vertices ^ r and £ eT is manifest. Effectively, as it is shown in table 
(|4*|). decays involving the £ eT vertex, are highly suppressed respect to the current instrumental sensitivity. 
The constraints on the rest of LFV vertices are much milder. However, current prospects to increase the 
statistics concerning LFV decays could improve such bounds significantly. 

In addition, we constrain some FC couplings in the quark sector by using experimental limits on leptonic 
B® decays as well as the B® — B d mixing. The leptonic B decays provide constraints on the products of 
lepton and quark FC couplings. On the other hand, by assuming that the charged Higgs boson is sufficiently 



-Fr<> 



heavy, the B® — B d mixing can be used to constrain the vertex ^b- We point out that although the B d 



b: 



mixing imposes severe restrictions to this vertex, the upper limits for £^6 could still give a sizeable and even 
dominant contribution to the decay B d —> fip respect to the SM contribution. Consequently, this vertex of 
the first generation can still be important for phenomenological calculations. 
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DINAIN. 
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